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Abstract 

On a laboratory-scale testing platform of impinging entrained-flow gasifier with four opposed burners, the flame images for diesel 
combustion and gasification process were measured with a single charge coupled device (CCD) camera. The two-dimensional 
multifractal detrended fluctuation analysis was employed to investigate the multifractal nature of the flame images. Sound power- 
law scaling in the annealed average of detrended fluctuations was unveiled when the order q > and the multifractal feature of 
flame images were confirmed. Further analyses identified two multifractal parameters, the minimum and maximum singularity 
Qmin and Qmax, scrviug as characteristic parameters of the multifractal flames. These two characteristic multifractal parameters 
vary with respect to different experimental conditions. 
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1. Introduction 

Coal gasification is an important chemical process 
for coal and the key to realize coal clean utilization. 
Entrained flow gasification is leading among coal 
gasification technologies and represents one of the 
cleanest ways of coal utilization. The entrained-flow 
gasification technology has been extensively applied 
to the production of ammonia, methanol, acetic acid, 
other chemicals and power generation through Inte- 
grated Gasification Combined Cycle (IGCC). The 
gasification process of an entrained-flow gasifier is 
very complicated, because it relates to the fluid flow 
under the condition of high temperature, high pres- 
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sure and heterogeneous state. Impinging stream flow 
configurations are characterized by streams of fluid 
jets impinging against each other in a confined vessel, 
which have proved useful in conducting a wide array 
of chemical engineering unit operations and enhanc- 
ing heat and mass transfe r between ph ases due to the 



high transfer coefficients (ITamin, 119941) . The opposing 



tensively studied both practically (Tamiret al. , 
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nat- 
The 



ural and social sciences ([Mandelbrot 
fractal behaviour of turbulent premixed fla me fronts 
has been manifested in man y experiments ( !Gouldii], 



19871: iGouldin et all Il989l: iMuravama and Takeno 
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1989 
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Wu et al.jl991l:lGoix and Shepherd,1993i:iYoshida et al 



1994allbl: LSmallwood et all Il995t lErard et all 1 1996 : 
Das and Evansl 1997 ). These studies provided a quan- 
titative description of flames, which enables us to better 
understand the dynamics of different types of com- 
bustion. On the other hand, to the best of our knowl- 
edge, the multifractal nature of impinging flames has 
not been studied. We shall adopt the two-dimensional 
multifractal detrended fluctuation an a lysis ( MFDFA) 
recently developed by Gu and Zhou ( 20061) for this 
purpose. The 2D MFDFA is a ge neralization of the 
ID DFA and MFDFA inven ted by IPeng et al.l (1 1994 
and iKantelhardt et al.l (l2002h . which is widely used to 
investigate fractal behavior of time series. The fam- 
ily of DFA approaches have the advantages of easy 
implementa tion, high precision, and low cornputa 
tional time jPeng et al. ,ln 
Gu and Zhoui l2006h . 



1994; Kantelhardt et al., 2002 



The paper is organized as follows. Section|2]reviews 
the algorithm of the two-dimensional MFDFA in detail. 
Under each experimental condition, we have recorded 
many flame images. Rather than analyze the images 
one by one under the same experimental condition, we 
propose to perform annealed averaging over dozens of 
images. The method of annealed averaging is also de- 
scribed. Section |3] outlines the schematic diagram of 
the experiment setup and the details of the experiments. 
Section |4] performs multifractal DFA analysis of the 
flame images, investigates the dependence of multifrac- 
tal characteristics with respect to experimental condi- 
tions, and discusses the physical interpretation of mul- 
tifractality. Section |5] summarizes. 

2. Methodology 

2.1. Two-dimensional MFDFA 

The ide a of DFA was invented originally by 



Peng et al.l (119941) to investigate the long-range de- 



pendence in coding and noncoding DNA nucleotide 



sequences and then generalized by IKantelhardt et al 
(l2002i) to study the multifractal nature hidden in 
time seri es, termed as multifr actal DFA (MFDFA). 
Recently, iGu and Zhou (l2006h generalized the one- 
dimensional MFDFA to higher-dimensional version, 
which is capable of analyzing multifracta l properties of 
higher-dimensional objects. According to lGu and Zhou 
(I2OO61) . the two-dimensional MFDFA consists of the 
following steps. 



Step 1: Consider a flame image t, which is de- 
noted by a two-dimensional array X{t\ i, j), where i = 
1,2,--- ,M, j = 1,2,-- - ,iV, and t = 1,2,--- ,T. 
The surface is partitioned into MgNs disjoint square 
segments of the same size s x s, where = [M/s] 
and Ng ~ [N/s]. Each segment can be denoted by 
Xy^w{t) such that ^ X{t]li+i, h+j) for 

^ ^ ii j ^ s, where li ~ {v — \)s and I2 ~ {w — l)s. 

Step 2: For each segment Xy^^{t), the cumulative 
sum Uy^w{t] i,j) is calculated as follows: 

i j 

Uv,wit;i,j) = X] 5Z ^v,wit;ki,k2) , (1) 

ki = l k2 = l 

where 1 ^ j ^ s. Note that u^^w itself is a surface. 

Step 3: The trend of the constructed surface Uu.u, (<) 
can be determined by fitting it with a prechosen bivari- 
ate polynomial function u. In this work, the following 
polynomial is adopted, 

Mi,,tu (t; «, i) = ai^ + bf + cij + di + ej + f , (2) 

where 1 ^ i,j ^ s, and a, b, c, d, e, and / are free 
parameters to be determined. These parameters can be 
estimated easily through simple matrix operations, de- 
rived from the least squares method. We can then obtain 
the residual matrix 

The detrended fluctuation function F{t; v, w, s) of the 
segment Xy,.w{t) is defined via the sample variance of 
the residual matrix as follows 

^ i=i j=i 

Step 4: The overall detrended fluctuation is calculated 
by averaging over all the segments, that is. 



Fg{t;s) 
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where q can take any real value except for g = 0. When 
g = 0, we have 



Foit;s) 



exp ■ 



v—l w—1 J 



(6) 

according to L'HopitaFs Rule. 

Step 5: Varying the value of s in the range from 



6 to Sr, 



min(M, N)/A, we can determine 



the scaling relation between the detrended fluctuation 
function Fq{t] s) and the size scale s, which reads 

. (7) 



2 



Usually, the scaling laws hold in a properly determined 
scaling range. 

Since AI and N are often not a multiple of the seg- 
ment size s, two orthogonal strips at the end of the 
profile may remain. In order to take these ending parts 
of the surface into consideration, the same partitioning 
procedure can be repeated starting from the other three 
corners. In this way, we have AMgNg segments and the 
overall detrended fluctuation is calculated by averaging 
over them. 

The final outcome of the MFDFA analysis is a fam- 
ily of scaling exponents h{q) which is a decreasing 
function of q for multifractal surface and remains con- 
stant for monofractals. I n the standard multifractal for- 
malism of lHalsev et al. I (ll986h based on partition func- 
tions, the multifractal nature is characterized by the 
mass exponent rig ) , whic h is a nonlinear function of q 



Kantelhardt et al.l (120021) showed that, for each q, we 



can obtain the corresponding traditional r((j) function 
through 



r{q) = qh{q) - Df 



(8) 



where Df is the fractal dimension of the geometric 
support of the multifractal measure. In this work, we 
have Df = 2. Res orting to the Legendre transform 
(iHalsev et al. . 119861) . We can also determine the singu- 
larity strength function a{q) and the multifractal spec- 
trum f{a) as follows 



a{q) = h(q) +qh'{q) 
f{a) = qa(q) - T{q) 



(9) 
(10) 



2.2. Annealed averaging 

To achieve better accuracy and higher statistical sig- 
nificance, we perform annealed averaging over 50 im- 
ages for each experimental condition, which are selected 
randomly from a huge data base. The annealed averag- 
ing gives the means of Fq{s) over the selected images. 
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(11) 



When (7 = 0, -Fo(s) can be calculated according to the 
following expression 

1 

ln[Fo(s)] = -^hi[Fo(t;s)] . (12) 



The characteristic functions of the ensemble of multi- 
fractal images can be determined according to Eqs. (|7} 



The shape and width of /(a) curve contain significant 
information about the singular flame images. In general, 
the spectrum has a concave downward curvature. We 
obtain two characteristic quantities, the width of mul- 
tifractal spectrum Aa = ckmax — amin and the differ- 
ence of fractal dimensions A/ = /(amax) — /(amin) 
of the minimum probability subset with a = c^max and 
the maximum subset with a = amin- These two quan- 
tities are used in this work to characterize different tur- 
bulent flames under different conditions. 



3. Experimental equipment and procedure 

The schematic drawing of the experimental apparatus 
is shown in Fig. [T] The maximum values of the opera- 
tion pressure and temperature were 1 MPa and 1500°C, 
respectively. The gasifier was cylindrical, vertically ori- 
ented. The combustion chamber was composed of a 15 
mm thick cast refractory shell and its inner diameter 
and length were 300 and 2200 mm, respectively. The 
cast refractory shell, wrapped with a 235 mm thick and 
low thermal conductivity fiber blanket to reduce the heat 
transfer, was protected by a stainless-steel column shell 
of 0.8 m in diameter and 2.5 m in height. Ports were 
located on the wall of the gasifier for viewing, tem- 
perature measurement, and insertion of a water-cooled 
probe. Opposed turbulent flow fields were obtained by 
four opposed round burners composed of inner and outer 
channels. Oxygen was fed into the outer channels of 
burners by steel cylinder, with a pressure-reducing valve 
to avoid pressure oscillations and achieve steady flow. 
The gas flow rates were measured by mass flow me- 
ters (D07-9C/ZM, Beijing Sevenstar Huachuang Elec- 
tronic Co., Ltd). Diesel oil was fed into the inner chan- 
nels of burners by a gear pump (A-73004-00 jj, Amer- 
ica Cole-Parmer Company), whose flux was determined 
gravimetrically with an elapsed timer and an electronic 
weight scale. 

In the gasification process, four burners were used to 
produce opposite jets of fuel that impinge at the center 
of the combustion chamber. The size of the burners is 
shown in Fig. |2] High relative velocities between the 
particulate matter and the gaseous phase in the central 
area provided good conditions for active diffusion and 
convection at the particle surface, and the high tem- 
perature together resulted in fast burning and gasifica- 
tion reaction under highly reducing conditions to pro- 
duce raw syngas. High-temperature gaskets interfaced 
the furnace segments and eliminated all leakage. From 
the reaction chamber, the raw syngas flowed into the 
quench chamber, where the raw syngas was cooled and 
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Fig. 1. Schematic diagram of experimental apparatus: 1— Liquid 
tank, 2— O2 steel cylinder, 3— N2 steel cylinder, 4— Pump, 5— Gas 
mass flow meter, 6— Burner, 7— Vortex flowmeter, 8— Slag discharge, 
9"Flame monitoring system. 

partially scrubbed by water and then discharged. The 
impinging flames are recorded with a flame monitor- 
ing system, fixed on the top of the gasifier. After each 
experiment, nitrogen was fed into inner channels by a 
steel cylinder to clean the burners. 



O2 Diesel O2 



Water 




Fig. 2. General view of the burner (all dimensions in mm) 

The schematic structure of flame image detector is 
shown in Fig. |3] The flame image detector is fixed on 
the top refractory wall of the gasifier, and the whole 
flames are visible during the experimental process. The 
flame monitoring system consists of a lens, an optical 
probe, a CCD camera (Panasonic WV-CP470), and a 
microcomputer. The CCD camera and its accessories 
are cooled by water to avoid overheating. The objec- 
tive lens is fixed at the front end of the optical probe 
and its surface is kept clean against dusts by a jet of 
nitrogen. The light conveyed by the optical probe en- 
ters the CCD camera. The camera has a 1/3-inch high 
resolution progressive scan interline-transfer CCD sen- 
sor with an array of 720 x 576 pixels and each pixel 



contains 3 bytes (one for each of the fundamental col- 
ors: red, green and blue). The video signal is transferred 
from the CCD camera to the microcomputer and stored. 

Power input Cooling Water-cooled Sweeping 

connector CCD camera water outlet jacket N2 pipe 



I 



Video out Purging 
connector nitrogen 



water inlet 



Optical 
probe 



Fig. 3. (Color online) Schematic structure of the flame image detector 

A typical flame image is illustrated in Fig. |4] The 
image analysis included the following steps. First, the 
video signal was digitized into 8-bit (256 gray values 
for each pixel) two-dimensional digital color images at 
a rate of 24 frames per second. Second, for each ex- 
perimental condition, 50 images were analyzed, each 
of them was converted into grey images with the level 
ranging from (black) to 255 (white). Third, in order 
to minimize the possible generated effects of the gasi- 
fier wall background, we cropped the grey images to 
remove most of the gasifier wall background but pre- 
served the whole flame resulting in 200-by-200 images. 
Four, multifractal analysis was then performed on the 
50 resultant images. 




Fig. 4. (Color online) The typical four-burner impinging flame image 



4. Results and discussion 

4. 1 . Annealed multifractal DFA of images 

For each experiment, we performed annealed mul- 
tifractal DFA on 50 arbitrarily chosen flame images. 
Fig. |5] shows the dependence of the annealed average of 
the detrended fluctuations In [Fg(s)] with respect to the 
scale s for six different orders q. The nice linearity of 
the lines indicates power-law scaling between Fq{s) and 



4 



s. The scaling ran ge spans abou t 2.15 orde rs of magni 



tude. A ccording to lMalcai et al.l (119971) and lAvnir et al 



([19981), the scaling range width for experimental frac- 
tality is about 0.5 to 2.0 orders of magnitude. It means 
that the power-law scaling observed in our experiments 
are quite sound. 



12 
10 

3 « 

^. 6 





— 9=0.6 




^9=1.4 




9=2.2 




—9=3.0 




^9=3.8 




^9=4.6 







1 2 3 4 5 

ln(s) 

Fig. 5. The annealed average of detrended fluctuations Fq{s) as a 
function of the scale s for six different values of q 

According to Eq. the slopes of the straight lines 
obtained by least squares regression of (s)] against 
ln(.s) in Fig.|5]give the estimates of the scaling exponent 
h{q). In Fig.|6]is illustrated h{q) as a function of q for 
< q ^ 5. We observe that h{q) is a nonlinear func- 
tion of the order q, which is the hallmark of multifrac- 
tality in the flame images. We note that the estimate of 
Fq{s) becomes statistically less significant for larger q 
due to the finite size of the flame images. It is also worth 
stressing that, no evidence of power-law behaviour is 
observed for q < 0. In other words, the scale invariance 
is destroyed for negative q. This phenomenon is not 
unusual in the multifractal analysis of experimental re- 
sults or natural measures. Examples include the growth 
probabilities at per imeter sites of diffusion -limited ag- 
gregations (DLA) (ILee and Stanleyl 119881) . the spatial 
distribution of the secondary electrons on solid surfaces 
and in the bulk tLi et al. 1995, 1 996i), the grow t h prob - 



ability of a solid-on-solid model dWang et al.L Il995h , 
the TEM imag es of four-layered GeAl film after laser 
irradiation (Sanchez et al.l Il992h . to list a few. 

The values of a{q) and /(a) can be computed nu- 
merically based on the Legendre transform ( 19110b . The 
resulting multifractal curve f{a) is plotted in Fig. |7] 
with respect to a. Two fundamental quantities amin and 
ttmax are determined, which characterize respectively 
the minimal and maximal singular sites of the turbu- 
lent flames. Moreover, the width of singularity strength 
Aa = Q!niax — Oinun Can be used as a measure of het- 
erogeneity of the singular flames. Similarly, /(ofmin), 
), and A/ = f{a 

max ) - ./(amin) can also be 




Fig. 6. Dependence of h{q) with respect to q 

used to quantitatively characterize the multifractal na- 
ture of the flames. 




Fig. 7. Multifractal spectrum for the flame images 

4.2. Relationship between multifractal parameters and 
experimental conditions 

In our experiments, we have investigated 95 con- 
ditions by varying the burner exit velocities of diesel 
and oxygen, M(m/s) and V{to? /s). Under each ex- 
perimental condition, 50 flame images have been used 
to calculate the ensemble multifractal parameters amin, 
Qfmax, Aa, /(a,nin), /(ckmax), and A/. Here we in- 
vestigate the possible dependence of these multifractal 
parameters with respect to experimental parameters M 
and V to identify characteristic multifractal parameters. 
Specifically, we adopt linear models as follows 



MY ^ Po + PiM + f32V 



(13) 



where MF stands for the six individual multifractal pa- 
rameters. We investigate three different types of models 
by posing /3i = or /32 = or freeing /?i and (32- This 
gives 18 models in total. For each model, we use the 
F-test to check if the model is statistically significant. 
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If the model is significant, we employ further the t-test 
to see whether the model coefficients are significantly 
different from zero or not. 

As a first step, we consider the six monovariate mod- 
els with (3i ~ 0. The F-tests show that the four models 
for Aa, /(amin), /(ctmax) and A/ are not significant, 
where all the p- values are greater than 13%. On the other 
hand, the p-values for a,„i„ and a,„ax are both 0.1%. 
This implies that both the multifractal parameters amin 
and ttniax are linearly dependent on the burner exit ve- 
locity V of oxygen. In addition, according to the t-tests, 
we find that the p-values of (3o and ^2 in these two mod- 
els are not greater than 0.1%. In other words, the coef- 
ficients (3q and P2 are significantly different from zero. 

Then, we investigate the six monovariate models with 
(32 — 0. The F-tests show that the five models for 
ttniin, Aa, /(omin), /(amax) and A/ are not signifi- 
cant, where all the p-values are greater than 29%, while 
the p-values for amax is 4.5%. This means that the mul- 
tifractal parameter Qmax is linearly dependent on the 
burner exit velocity M of diesel, significant at the 5% 
level. In addition, according to the t-test, we find that 
the p-values of Pq and /3i in this model are not greater 
than 5%. In other words, the coefficients Pq and Pi are 
significantly different from zero. 

Now we turn to study the six bivariate models (l3[ 
with different dependent variables MF. The resluts are 
listed in Table[T| According to the F-tests, the four mod- 
els for Aa, /(amin), /(amax) and A/ are not signifi- 
cant whose p-values are greater than 14%. For amin and 
ttmax, the p-values obtained from the F-tests are less 
than 1%. Therefore, both the multifractal parameters 
amin and amax are linearly dependent on the velocities 
of oxygen and diesel. However, the t-tests shows that 
the coefficients Pi for both models are not significantly 
different from zero. Speaking differently, the multifrac- 
tal parameters amin and amax depend strong upon the 
burner exit velocity of oxygen and weakly on the ve- 
locity of diesel. Furthermore, the minimal and maximal 
singularity strengthes amin and amax are characteristic 
multifractal parameters that can be used to quantify the 
multifractal nature of the turbulence flames. 



4.3. Discussion 

We have confirmed that the impinging flame images 
exhibit multifractal properties. A natural question arises 
asking which processes lead to such multifractality. Sev- 
eral examples taken from other fields in literature pro- 
vide c lues for addressing this question. iGodano et al 
(Il997h argued that the multifractal nature of the tempo- 



Table 1 

Identification of cliaracteristic multifractal parameters using linear 
model dD 





Coefficients 


t-test 


F-test 


MF 


Po Hi A 


PO Pi P2 


F p 


^min 


2.4734 0.4538 -0.0027 


0.000 0.589 0.001 


5.576 0.005 




3.1728 -0.6503 -0.0013 


0.000 0.117 0.002 


7.477 0.001 


/(^min 


) 0.3067 2.1950 -0.0011 


0.649 0.260 0.555 


0.729 0.485 




) 1.9911 -0.0576 0.0002 


0.000 0.610 0.118 


1.279 0.283 


Aa 


0.6992-1.1035 0.0014 


0.023 0.208 0.094 


1.942 0.149 


A/ 


1.6845 -2.2529 0.0013 


0.013 0.243 0.491 


0.822 0.443 



ral clustering of earthquakes can be interpreted in terms 
of diffusive processes of stress in the Earth's crust. In 
order to explain the multifractal p ropert ies in rainfall 
data, lOlsson and Niemczvnowicz (Il996h made an as- 
sumption that a large-scale flux is successively broken 
into smaller and smaller ones in cascades, each receiv- 
ing an amou nt of t he tota l flu x specified by a rn ultiplica- 
tive process. iLed ( 120021) and lLee et alj ( l2003h also em- 
ployed a random multiplicative process to understand 
the multifractal characteristics in air pollutant concen- 
tration time series. In the present case, we submit that 
a possible explanation for the physical origin of mul- 
tifractality in flame images is the atomization and im- 
pinging processes. 

Due to the complexity of the atomization process, it 
is difficult to clearly describe the mechanism and also 
impossible to combine all the influencing factors into 
one model, such as the equipment dimensions, the size 
and geometry of burner, the physical properties of the 
dispersed phase and th e continuou s phase , and the op- 
erating mode. Indeed, IZhou et al. proposed a 
stochastic multiplicative cascade model for drop break- 
down in atomization processes, which works we l l in th e 



prediction of drop size distribu tion dLiu et al. . l2006h 



Moreover, IZhou and Yul (1200 lb studied the multifrac- 
tal nature of drop breakup in the air-blast burner atom- 
ization process. They applied the multiplier method to 
extract the negative and the positive parts of the /(a) 
curve with the data of drop-size distribution measured 
using dual particle dynamic analyzer They proposed a 
random multifractal model with the multiplier triangu- 
larly distributed to characterize the breakup of drops, 
the agreement of the left part {q > 0) of the multifrac- 
tal spectrum between the experimental result and the 
model is remarkable. Hence, the breakdown of diesel 
drops in the gasifier follows a multifractal process. 
Four equal suspension streams flow against one an- 
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other at high velocity (> 35 m/s) and impinge at the cen- 
ter of the gasifier, resulting in a highly turbulent zone. 
The gas flows decrease their axial velocity down to zero 
at the impingement plane and then disperse radially, 
while particles penetrate back and forth between the op- 
posed streams driven by inertia and friction forces. This 
impinging process leads to locally singular distributions 
of the raw materials. 

5. Conclusion 

The impinging process has proved to be one of the 
most effective methods enhancing heat and mass trans- 
fer in multiphase environment. On a laboratory-scale 
testing platform of impinging entrained-flow gasifier 
with four opposed burners, the flame images for diesel 
combustion and gasification process were measured 
with a single charge coupled device (CCD) camera. 
Ninety-five experimental conditions were investigated. 

The multifractal properties of the turbulent flames 
have been investigated using the two-dimensional multi- 
fractal detrended fluctuation analysis, which is accurate 
and easy to implement for image analysis. Nice power- 
law scaling is unveiled in the annealed average of de- 
trended fluctuations when the order q > 0. The scaling 
exponent h{q) is a nonlinear function of q, which con- 
firms the multifractal nature of the flames under investi- 
gation. We argue that the multifractality in gasification 
flames stems from the multiplicative process of atom- 
ization and the impinging process of raw materials. 

We analyzed the relationship between six multifrac- 
tal parameters (amin, "max, Aa, /(amin), ./("max) and 
A/) and the velocities of oxygen and diesel {V and M). 
Two multifractal parameters amin and Ofmax have been 
identified by extensive F-tests as characteristics of the 
observed multifractal nature, which are dependent on 
the velocities linearly. The t-tests show that the velocity 
of oxygen has greater impact on the multifractahty of 
flames. These analyses enable us to gain a better under- 
standing of the complexity of the combustion dynamics 
of four-burner impinging entrained-flow gasification. 
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Notations 

s scale of boxes 



T = 50 


numbGr of images 


xa i) 


two-dimensiona.1 array 


"'VjW V'') ^1 J ) 


cumulative sum for image t 


ti „, „,,('} ')] 
'-''V \ J ! 


fitting bivariate polynomial 


'") J ) 


residual matrix 


n 
H 


order of detrended fluctuation function 


Fr,(s\ 


detrended fluctuation function 




fractal dimension 


h(a) 


scaling exponent function 


T(a) 


mass exponent function 


a(q) 


QiTiCTn 1 nn ctrpntrth fiinr'tinn 


ttmax 


maximum singularity 




minimum singularity 


Aa 


width of multifractal spectrum 


/(«) 


multifractal singularity spectrum 


A/ 


difference. A/ = /(omax) - /(amin 


M 


burner exit velocity of diesel (m/s) 


V 


burner exit velocity of oxygen (m/s) 


/9o,/3i,/32 


model coefficients 


P0,Pl,P2 


p-values of model coefficients 
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